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We have studied scalar perturbations as well as fcrmion perturbations in pure de Sitter 
^ ■ space-times. For scalar perturbations we have shown that well-defined quasinormal modes 



in d-dimensions can exist provided that the mass of scalar field m > The quasinormal 
modes of fcrmion perturbations in three and four dimensional cases have also been investi- 
gated. We found that different from other dimensional cases, in the three dimensional pure 

^ ■ de Sitter spacetime there is no quasinormal mode for the s-wave. This interesting difference 

!>■ 

' caused by the spacial dimensions is true for both scalar and fcrmion perturbations. 

o ■ 
o 
o 
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INTRODUCTION 



It is well known that the surrounding geometry of a black hole will experience damped oscilla- 



^ . tions under perturbations. The frequencies and damping times of the oscillations are entirely fixed 

\ by the black hole parameters and independent of the initial perturbations. These oscillations are 

called "quasinormal modes" (QNM), which is believed as a characteristic "sound" of black holes 
and would lead to the direct identification of the black hole existence through gravitational wave 
observation to be realized in the near futureO]. Due to the potential astrophysical interest, a great 
deal of effort has been devoted to the study of black holes' QNMs. Most of these studies were 
concerned with black holes immersed in an asymptotically flat spacetime P]. Considering the case 
when the black hole is immersed in an expanding universe, QNMs of black holes in de Sitter(dS) 
space have also attracted much attention 
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Motivated by the discovery of the AdS/CFT correspondence, the investigation of QNM in anti- 
de Sitter(AdS) spacetimes became appeahng in the past several years. It was argued that the QNMs 
of AdS black holes have direct interpretation in term of the dual conformal field theory(CFT) 4|- 
In dS space the relation between bulk dS spacetime and the corresponding CFT at the past 
boundary and future boundary X+ in the framework of scalar perturbation spectrums has also 
been discussed A quantitative support of the dS/CFT correspondence was provided. 

Recently QNMs in asymptotically flat spaces have acquired further attention, since the possible 
connection between the classical vibrations of a black hole spacetime and various quantum aspects 
was proposed by relating the real part of the QNM frequencies to the Barbero-Immirzi(BI) pa- 
rameter, a factor introduced by hand in order that loop quantum gravity reproduces correctly the 
black hole entropy ll| . In order to see whether this quantum connection is true in Schwarzschild 



dS (SdS) spacetime, a number of extentions have been made|l2|| [14l| . For the nearly ex- 
treme SdS black hole, like the Schwarzschild spacetime, the real part is found still proportional to 
the black hole surface gravity, but instead of the integer n labelling the modes, with proportional 



coefficients \/^(^ + 1) ~ i foi' scalar and electromagnetic perturbations and J + 2)(^ — 1) — 4 



for gravitational perturbations. It is too ear 
^-dependent real part to the BI parameter 



y to confirm the recent conjecture by relating the 
12l |. The imaginary part of the QNM frequencies for 
the nearly extreme SdS black hole was found having an equally spacial structure with the level 
spacing depending on the surface gravity of the black hole, which agrees to the result in asymptot- 
ically flat spacetime and is independent of whether the perturbation is scalar, electromagnetic or 
gravitational 



12l |. This result was confirmed even for the very small SdS black holes |l^ and later 
further supported by using the Born approximation [l4 1 . 

The dependence of the surface gravity for both the real part and imaginary part of QNM fre- 
quencies suggests that there is a possible connection between the QNMs and thermodynamics of 
the black hole horizon. But what is the effect of the cosmological horizon here? It is believed that 
cosmological horizon has very similar thermodynamical behavior to that of the black hole horizon 
[l^ . Is there connection between the QNMs and surface gravity of the cosmological horizon? In 
order to answer this question, in this paper we are going to investigate the pure dS spacetime. It 
was argued that for a massless minimally coupled scalar field, there exists no QNMs in the pure 
dS spacetimes, however for a massive scalar field, there do exist well-defined QNMs^^. We will 
examine this arguement in different topological pure dS spacetimes by considering scalar perturba- 
tion and fermion perturbation. Different from the purpose to test the dS/CFT correspondence^^, 
we here concentrate on the region within the cosmological horizon. 
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The structure of the paper is as follows: In Sec. II, we will study the scalar perturbation in 
d-dimensional pure dS spacetime with a spherically symmetric space. The extension to topological 
pure dS spacetimes will be shown in Sec. III. In Sec. IV, we will present the discussion of the fermion 
perturbation. Our main results will be summarized in Sec.V. 



II. SCALAR PERTURBATION IN PURE DE SITTER SPACE 

The static metric of a d dimensional dS space reads 

ds^ = -f{r)dt^ + rHr)dr^ + r^dnl_^ , (1) 

where /(r) = 1 — r'^/l'^{l is the minimal radius of dS space) and r'^dQ,'^_2 represents the metric on 
the d — 2 dimensional sphere S'^"^ of radius r. 

We begin our discussion with a massive scalar field <I>, satisfying the Klein-Gordon equation 

^■'".^ = (2) 

This equation can be separated by $ = (u(r)/r~2~)e^"^*Y£(r2^_2)- Here the spherical harmonic 
Yi{Q(i-2) is the eigen function of d — 2 dimensional Laplace-Beltrami operator V^_2 with the 
eigenvalue + d — 3). Using the tortoise coordinate, = J dr//(r) = / tanh^^(r//), we can 
write the radial part into a Schrodinger-like equation 

^ ^ + V{r^)u = u?u (3) 



with effective potential 



A B 

where A and B are defined by 

d-2.d-2 . 2,2 
A = ^(^ + l)-m¥, 

B = i{i + d-3) + ^{^-l). 

When B > 0, the effective potential @ diverges to a positive infinity at the origin(r = 0) and 
vanishes at the cosmological horizon. On the other hand, B < 0, the potential falls down to 
negative infinity at the origin which indicates the instability of the perturbation. We will show it 
in detail later. 
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In the form of a new variable 2 = 1/ cosh^(r*//), equation 1^ becomes 



z{l - z)u" + {1 - -z)u' + - 



,2;2 



B 



1-z 



+ A 



u = 0. 



(5) 



Further using the ansatz u = — z)^F{z), we have 



z{l-z)F"+ 1 + 2q- (2a + 2/3 + -)z 
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F = 



(6) 



If we properly select values of a, (3 to make terms like ^ and disappear, the solution to equa- 
tion® is the standard hypergeometric function. It is not hard to find that ~ exp(ib iojr^) when 
z 0(that is, approaching the horizon), so that the two independent solutions exactly correspond 
to the incoming and outgoing waves at the cosmological horizon. The general solution is 

U{z) = Ul{z) + U2{z) 

= Ci2""(l - z)^ 2i^i(a - c + 1, 6 - c + 1, 2 - c, z) 
+ C2z"(l-z)^ 2Fi{a,b,c,z) 

where Ci, C2 are constants, and the parameters of the hypergeometric function are given by 



(7) 



2q + 1 



a + P + -{I + 



b = a + P + -{I - VTT4A) 

Note that the condition of writing the solution in the form of (7) is that c is not an integer. 

We will restrict ourselves in the well-accepted definition of the quasinormal mode which is, 
in this special case, that these modes are purely outgoing waves at the cosmological horizon and 
vanish at r = 0. Indeed this boundary condition is determined by the behavior of the effective 
potential. Since the sign of a is arbitrary, it is easy to check that choosing either ui or U2 as 
the outgoing wave does't lead to any difference. Here we take Ci = 0(so a = — za;//2)jhence the 
incoming wave is eliminated. According to the property of hypergeometric function f2g|. we can 
change the wave function into 



U2iz) = C2Z-il-zf^^^^^^^^^^^ 



2Fi{a, b,a + b - c + 1,1 - z) 



r(c-a)r(c-6) 

2-Fi(c — a, c — b,c — a — b + 1,1 



T{a)m 



(8) 
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To make it zero at the origin, we use the poles of Gamma functions. And to obtain the discrete 
poles which result in the level spacing frequencies we should assume that 

1 B 

P{^-P) = -^<o, 

which is identical to our previous analysis of the effective potential. It is interesting to find that 
in three dimensional case, this condition is equivalent to P > 1/4, which implies that there exists 
no QNM for s-wave perturbation. 

There are two sets of poles, (i) When /3 = ^^^^ — ^ > ^ (note that for three dimensional 
case we have excluded the s-wave perturbation for the reason mentioned), the poles are o = —n 
or b = —n (n = 0, 1,---). These poles indeed remove the divergent part of the wave function 
around the origin since both c and a + b — c = could not be nonpositive integers and the 

numerator remains regular. (ii)When (3 = — + 1 < 0, the poles are c — a = —n or c — 6 = —n 
(n = 0, 1, • • • ). Because both c and c — a — b = could not be nonpositive integers, these 

poles are well defined. 

The corresponding frequencies are 

io = -j{2n + £ + h±) (9) 
or uj = -j{2n-£-d + 3 + h±) (10) 

where /i± = ^ - ^{^Y-mH'^. 

The above discussion should be restricted to the condition that c is not an integer. Now let us 
see what will happen when this condition is violated. It only occurs in the massless case m = 0. 
We assume that c = k, where k = 0, ±1, ±2, • • • . When c = A; > 1, we can obtain the solution 
satisfying the boundary condition 



_iuA . , d+2e-2 

u{z) oc z 2 [1 — z) 4 2Fi{a,b,k, z) (11) 



and the corresponding quasinormal frequencies are 



a; = ^(fc-l) (12) 



When c = k <0, the proper solution is then 



, . iuil , . d+2t—2 , ^ , ^ 

u{z) oc Z-— {I- z)^^ 2Fi{a-k + l,b -k + 1,2- k,z) (13) 



and the quasinormal frequencies are 



u: = -'-{k-l). (14) 
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It is surprising that all these modes © (ITU]) (IT^ ( (11)) are just some kinds of "distribution" along 
the radius. To show that, we compute its flux(|]^ 

T = ^\^XR*drR - RdrR*) (15) 

where R = U2{r)/r^~ is the radial factor of the wave function. If frequencies w in ® (|1U()(|12I)(|14|) 
are all purely imaginary values, a = —iujl/2 is real and so are the parameters a, b and c. Then 
R{r) is proportional to a purely real function defined in the region < r < I, which results in the 
vanishing of the flux ()15j) everywhere as well as at the horizon! This is in contradiction with the 
definition of QNM. 

However for the massive case, one can see from the scalar perturbation spectrum ((m) (|lUj) that if 
m'^P > (^^)^, the frequencies are not purely imaginary values which insure that the flux does not 
vanish. These modes are purely outgoing waves at the cosmological horizon and so are well-defined 
QNMs. Therefore there exists the lowest bound (m > (d — l)/2l) of the mass of scalar field that 
permits QNMs to survive. We rewrite the corresponding QNM frequencies as follows: 

^ = ±y [m'l' - i^f] ^ - ]i2n + ^ + ^) (16) 
or a; = ±i [mH' - {^f] ^ - ^(2n - £ + 3 - ^) (17) 

This result confirms the argument that QNMs cannot exist for massless scalar field Here 
we further present the lowest mass bound for the scalar field to posses QNMs. 



III. SCALAR PERTURBATION IN TOPOLOGICAL DS SPACE 

In AdS spacetime, the QNMs of different topological black holes have been studied recently. It 
was found that black hole topology influence a lot on the QNMs of scalar perturbations P|. Here 
we would like to extend the investigation to the topological dS space. 

The metric of the topological dS space is showed as with lq| 

„, , , 2Gm 

/(0 = ^+^-7J (18) 

where k = 0,-1 correspond to two kinds of hypersurfaces ^^^-2: Ricci fiat space R'^-^ and negative 
constant curvature space H'^~^ , respectively. Note that d > 4 since d = 3 is trivial. 

The wave equation for the scalar perturbation still has the form as but the effective potential 



Kn ^d-2fr d-2 d-2 f 2 
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where Kq is the eigenvalue of Laplace-Beltrami operator V^_2- It is easy to find that the effective 
potential tends to zero at cosmological horizon, while it goes to negative infinity when r — > 0. Due 
to the existence of negative effective potential for the scalar field equation of motion, the bound 

n 

states can be formed leading to the growing modes instead of decaying modes J,9|. 

The growing modes imply that these topological dS space are not stable. This result is not so 
surprising, since both of these topological dS spaces possess naked singularity. 

In AdS case, though the wave amplification behavior exists in hyperbolic space, for flat space it 
still exhibits decaying wave for scalar perturbation Q]. The growing modes for the flat hypersurface 
in dS case can be considered as one more difference between dS and AdS spacetimes. 



IV. FERMION PERTURBATION 
A. Three dimensional case 



We are interested 



The fermion perturbation in a BTZ background was studied in [3, 22, [2^. 
to generalize the investigation of the fermion perturbation in a pure dS space. There is a little 
difference between three and four dimensional case due to the different representations of gamma 
matrixes. We first consider the three dimensional case. After the coordinate transformation t = t, 
r = I cos fi and (j) = ip/l, metric (0) becomes 

ds^ = - sin^ n dr^ + /^d//^ + cos^ dcp"^. (20) 

We begin with the Dirac equation 

l^ea^id^ + r^)^' + m^ = 0, (21) 

where 7" are the conventional spin matrixes. In three dimensional case we choose the representation 
of them in term of Pauli matrixes: 7^ = ia"^, 7^ = cr^, 7^ = a^. The triads e°i/ are 

e''T- = sin^, e^fj_ = l, e^^ = cos (22) 

The spin connection, defined by T^, = ^'~f^"'^^^ea^eb\-i^, has only two non-vanishing components 

r. = -^cos^7'V, = -^sinMT'V'. (23) 

We write out equation (j^Q explicitly 

1 1 cosu s'mn . ia'^ r. r^ 1 „ 

7^ (^M + T— ^ - ^^-^) + - — + d^ + m^ = 0. 24 

/ 2 sm n 2 cos sm cos /i J 
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Separating the equation by 



we arrive at the radial part equation 



-iu)T ^—ii ip/l I \J/-|^ 



{d^ + -. ^1 = -{ml + *2, 

sm ^ cos ji 

{d^-——)^2 = -{ml - 



sm ^ cos ^ 

Introducing a new set of wave functions V'i)'02! which relate to ^'i, ^'2 by 



*i + ^'2 = (1+tanV) 4y/1 + aan// (V'l + V'2) 
*i - *2 = (1 +tan^/x)~^v^l - ztan// {ipi - ■(/'2) 



(25) 
(26) 



(27) 
(28) 



and using a new variable y = tan /i for the sake of convenience, we turn Dirac equations into 



:i + ?/Wi + (-^y + — 

y 

+ )V'2 

2/ 



(m/ + 2) + (^^ ~ 



(1 + 2/^)5j,V'2 - (- 



(ml + 



(i£ - Zw) 



^2, 



(29) 
(30) 



Further by setting z = —y^ and choosing the ansatz ipi = z"'{l — z)^F{z), from the coupled 
equation ()29()H3U() . we find the purely outgoing solutions in terms of hyper geometric function: 



^P^=B,z''{l-z)'' 2Fi{a,b, c, z), 



^2 = B2^2+"(l - z 

ah 



1/3 



a + (i) 2-^1 (a, 6, c, z) 



+— (1 - z) 2Fi{a + 1, 6 + 1, c + 1, z) 



where the constants a,/3 ,the hyper geometric parameters a,b,c and coefficients .61,-62 are 



(31) 



(32) 



iujl 

,1 im/. 



/5 = ±(4--), 
a = -- + a + p, 



1 + . 



+ a + /?, 



(33) 



5i = B2{--a±P). 
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We now consider the flux along the radial direction 



Using the property of hypergeometric function 
when z — oo 



i^2r). (34) 

we can obtain the asympotic behavior of ^\ 



^i = i?i(i-i)^(-ir{ 



z 

--r^'- ""rf^Kjl ,F,{a,a-c + l,a-b+l, k 
r(c — a)T[b) z 



~i?i(-ir 
Similarly for ^2 we have 

*-2~Bj(-l)°+4{ 



r(c — b)T[a) z ) 



T{c-a)T{b) 



T{c)T{a-b) 
T{c-b)T{a) 



r(c)r(a - b) 



4 " + r(c-a)r(6)+^ r(c-6)r(a) 



+- 



ab 



c L 



(-^) 



1 i.a+/3-a r(c + i)r(b-a) ,_ a+„+^_;, r(c + i)r(a - fe) 

r(c-a)r(6 + l)^ ' T{c-b)T{a + l) 



To make the flux T vanish at the origin, we should have poles of gamma functions 



c — a = — n, or 6+1 
c — b= —n, or a + 1 



-n, 



-n, 



for ^ > 
for i<0 



Thus the quasinormal frequencies are 



i 3 
LU = —m — -{2n + 



i 1 /IN 

or uj = m — + 2 "I" 

i 1 

or a; = m— -(2n i), 

I 2 



for ^ > 
for £ < 



(35) 



(36) 

(37) 
(38) 

(39) 
(40) 



r 2 

Prom H34 |) - (|36|) . we know that when i = the flux would not vanish at the origin, in contra- 
dictory to the QNM's definition. So it is true for the fermion perturbation as well as the scalar 
perturbation that there is no QNMs for the s-wave in the three dimensional pure dS space. 



B. Four dimensional case 



We now turn our discussion to the four dimensional case. The metric and the Dirac equation 
are showed as and H21|l but with d = 4. In this spherical coordinate, a general formalism was 
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provided in [20| (see also 



2J|). We will follow this setup. With the separation 



-iujt 



<f7 



where are the two component spinors with j = j^l) we arrive at the radial part equation 
. If we use a new variable r = I sin fj,, the equations are like 



id, + ^)G^ 



sm fj. 



-{ml 



ujI 



)F^ 



cos fl 



(41) 
(42) 



sm fj. cos fi 

where At+ and k_ are positive and negative integers. These equations resemble equations H25|)(|26|l 
in the form, if one replaces zca;/,£, ^'i, ^2 there with K±,iujl,G^, here. So the general solutions 
are similar. But things are a bit complicate here. We cannot find the purely outgoing solutions as 
we did in the three dimensional case. Thus we first consider the flux along the radial direction at 
the origin r = 0(or equivalently z = here) 



sine (IV^IV 



±|2|.^± |2 
jm\ 



(43) 



We can find the solutions that make the flux vanish at the origin. For the case (G"^, F"^, k+), they 
are 



^P+ = B+ z^{l - zf 2Fi{aiM,ci,z), 



B+ z2(\-zf 



1 , 1 + 
- + - ( 

LV2 + ^ 2 



iijjl 



+ /3)2:) 2^^1(01,^1, Ci,2:) 



+^2;(1 - z) 2Fi(ai + 1,61 + l,ci + \,z) 

Cl 



(44) 



(45) 



where the hypergeometric parameters a, h and c and the relationship between the coefficients 
are 

+ 1 — iw/ ^ , 2 + Kj_ + iujl 3 
oi = ^ Vli, bi = ^13, ci=K+ + -, 

For the case {G^ , F~ , k^), they are 



z 2 (1 - z)^ 2-^1(02, ^2, C2,Z) 



V^2- = ^2^ ^^(l_^)/3^(_i±^ + !|i_/3) 2Fi(a2,62,C2,z) 

+ 2Fi(a2 + l,62 + l,C2 + l,z) 

C2 



(46) 



(47) 
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where the hypergeometric parameters and coefficients are 

K_ + iujl 1 — K_ + iu)l 1 

02 = ^ h/3, &2 = ^ C2 = -K_ + -, 

iu)l 

B^ = B,i— + — ±f3). 

Note that /5 is the solution of equation + + = 0. 

To obtain the QNMs, we require that the waves are purely outgoing at the horizon. In other 
words, the wave function is of the form (— around the horizon (z — oo). To eliminate the 
ingoing part in the wave function ^^1.2! again, we appeal to the property of hypergeometric function 
illustrated in (j^SJ). The incoming and outgoing parts are then detached which enable us to remove 
the incoming one thoroughly by setting poles, which are 

ci - 6i = -n, or ai + 1 = -n, for (48) 

C2 — b2 = —n, or 02 + 1 = — n, for (49) 

The corresponding frequencies of QNMs are obtained as: 
For 

i 3 
a; = -m - y(2n + + -), or 

i 1 

uj = m- j{2n + K+ + -). (50) 



For {G-,F-,K- 



i 1 
(jj = —m — y(2ra — «:_ + -), or 



i 1 

u) = m — y(2ra — K- — -). (51) 

These results are very similar to that of three dimensional case, showing consistent behavior of 
the fermion perturbation in pure dS space. 

V. CONCLUSIONS AND DISCUSSIONS 

We have studied the QNMs of the scalar and fermion perturbations in the pure dS space. For the 
scalar perturbations, we have confirmed the argument(see ^|) that no QNM exists for the massless 
scalar perturbations in four dimensional space or, more generally , in arbitrary dimensional case. 
To allow the existence of QNMs of scalar perturbation, we have found that there is a constraint on 
the mass of scalar field, that is m > Moreover, we have found that in the three dimensional 
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pure dS space QNM does not exist for the s-wave scalar perturbation. This is a special result in 
three dimensional case and does not exist in other dimensions. 

We have also extended the discussion of scalar perturbation to the topological dS spaces. In 
two special cases, the flat hypersurface and the hyperbolic hypersurface, perturbations experience 
amplification behaviors. Here again we see that topology influence the behavior of the perturbation 
as that we observed in AdS situation 8]. In AdS cases, the wave amplification was only found in 
the hyperbolic spaces Q], however in dS spaces the growing modes are obtained in all topological 
spaces. This serves as an additional difference between the dS and AdS spacetimes. 

We have also investigated the fermion perturbations in the three and four dimensional pure dS 
spaces. The well-defined QNM frequencies were obtained. Similar to the scalar perturbation, we 
again found that for three dimensional case no QNMs exists for the s-wave perturbation. 

Examining the QNM frequencies of both scalar and Fermion perturbations, we found that 
though the real part of scalar perturbation is proportional to the surface gravity of the cosmological 
horizon, this dependence disappears in Fermion perturbation. In additional to the difficulty to fully 



13], it is also too early to 



understand the quantum connection between QNM and loop gravity 
see the relation between the QNMs and thermodynamics of horizons. 
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